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A combined analytical-computational method is developed to calculate the 
transient pressure field and dynamic coefficients for high-pressure annular seal 
configurations which may be used in interstage and neck-ring seals of multistage 
centrifugal pumps. The solution procedure applies to constant-clearance or 
convergent-tapered geometries which may have different (but directionally-homogeneous) 
surface-roughness treatments on the stator or rotor seal elements. It applies in 
particular to so-called "damper-seals" which employ smooth rotors and deliberately- 
roughened stator elements to enhance rotor stability. 

Hirs f turbulent lubrication equations are modified slightly to account for 
different surface-roughness conditions on the rotor and stator. A perturbation 
analysis is employed in the eccentricity ratio to develop zeroth and first order 
perturbation equations. The zeroth-order equations define both the leakage and 
the development of circumferential flow due to shear forces at the rotor and 
stator surfaces. The first-order equations define perturbations in the pressure 
and axial and circumferential velocity fields due to small relative motion 
between the seal rotor and stator. The solution applies for small motion about 
a centered position and does not employ linearization with respect to either the 
taper angle or the degree of swirl, i.e., the difference between the circumferential 
velocity at the given axial position and the asymptotic circumferential-velocity 
solution. 

Test results for four different surface-roughness confirm the predicted net 
damping increase for "damper seals". A round-hole-pattern stator yielded the 
highest net damping and lowest leakage of all seals tested. The seals are sub- 
stantially stiffer than predicted, but the theory does an adequate job of predicting 
net damping. 


NOMENCLATURE 


Dimensionless coefficients defined in Appendix A 

cT, (T Dimensionless damping coefficients defined by Eq. (34) 

f (z) Dimensionless clearance function defined by Eq. (9) 

The results reported herein were partially supported by NASA Contract NAS8-33716 
from the George C, Marshall Space Flight Center; Contract Monitor Frank Garcia. 
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h(z) = H/G Dimensionless clearance function 

h- First-order perturbation clearance function defined by 

Eqs. (11) and (18) 

k, K Dimensionless seal stiffness coefficients defined by Eq. (34) 

m, if Dimensionless mass coefficients defined by Eq. (34) 

ms, ns Dimensionless empirical turbulence coefficients for 

mr, nr stator and rotor 

p Fluid pressure (F/L 2 ) 

p* = p/pV 2 Dimensionless fluid-pressure introduced in Eq. (7) 

^ , p' Dimensionless fluid-pressure perturbations introduced in 

0 1 Eq. (11) 

q Taper-angle parameter defined in Eq. (XQ) 

t Independent variable time (T) 

\x = U„/V Dimensionless axial and circumferential velocity components 
introduced in Eq. (7) 

Uq = U q /R03 

u eo , u 01 a Zeroth and first-order perturbations in u 0 

U Z0 J U Z1 Zeroth and first-order perturbations in u^ 

z = Z/L Dimensionless axial coordinate 

A Test orbit amplitude (L) 

Dimensionless coefficients defined in Appendix A 

C Nominal seal radial clearance, CL) x 

C d Seal discharge coefficients defined by Eq. (16) 

C 0 , Entrance and exit clearances, respectively, (L) 

H(z, 0, t) Clearance function, illustrated in figure 2, and defined 
in Eq. (17), (L) 

H (z) Centered-clearance function defined by Eq, (9), (L) 

L Seal length (L) 

F g Seal supply pressure (F/L 2 ) 

AP Nominal pressure-drop across seal (F/L 2 ) 
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R 


Seal radius (L) 


R = 2pVH/]i Axial Reynolds number 
a 

R = 2pVC/|i Centered-position, axial Reynolds number 
ao 

T =i L/V Transit time for a fluid element to traverse the seal 

U^, Ug Axial and tangential bulk-flow fluid velocity components (L/T) 

V(z) Centered-position axial fluid velocity (L/T) 

X, Y Radial seal displacements (L) 

Z, R0 Spatial coordinates illustrated in Figure 2 

a Seal taper angle illustrated in Figure 3 

£ = e/g Seal eccentricity ratio introduced in Eq, (11) 

e = e/2- Relative roughness 

Li 

5 Inlet pressure-loss coefficient 

A , A Dimensionless stator and rotor friction-factors defined in Eq, (15) 

s j. 

T = t/T Dimensionless time 

oa Shaft angular velocity (T l ) 

Q Shaft precessional velocity (T *), introduced in Eq, (22)- 

INTRODUCTION 

Figure 1 illustrates the two seal types which have the potential for developing 
significant rotor forces. The neck or wear-ring seals are provided to reduce the 
leakage flow back along the front surface of the impeller face, while the interstage 
seal reduces the leakage from an impeller inlet back along the shaft to the back- 
side of the preceding impeller. Pump seals may be geometrically similar to plain 
journal bearings, but typically have clearance to radius ratios on the order of 
0.005 as compared to 0.001 for bearings. Because of the clearances, and normally- 
experienced pressure differentials, fully-developed turbulent flow normally exists 
in pump seals. 

As related to rotordynamics , analysis of seals has the objective of defining 
the reaction force acting on a rotor as a consequence of shaft motion. For small 
motion about a centered position, the relation between the reaction-force components 
and shaft motion may be expressed by 
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The off-diagonal coefficients in Eq. (1) are referred to as f, cross-coupled ,! and 
arise due to fluid rotation within the seal. Seals, unlike plain journal bearings, 
develop significant direct stiffness values K in the centered, zero-eccentricity 
position due to the distribution between (a) inlet losses, and (b) the axial 
pressure gradient due to wall-friction losses. Lomakin [1] initially pointed out 
the phenomenon. Both analysis [2] and experiments [3] have shown the Eq. (1) holds 
for fairly large eccentricies on the order of 0.5; i.e., the dynamic coefficients 
tend to be relatively insensitive to changes in the static-eccentricity ratio. 

Prior analytical and experimental developments have generally examined 
"smooth" seals where both stator and rotor elements of the seal are assumed to 
have the same nominally smooth surfaces . A- review of the analytical and experimental 
developments for this type of seal is provided in references [4] and [5] and will 
not be repeated here. The subject of this investigation is the so-called "damper- 
seal" configuration recently proposed by von Pragenau [6], which employs a smooth 
rotor and a deliberately surface-roughened stator element. For the same surface 
roughness on the rotor and stator, the asymptotic, circumferential, bulk-flow 
velocity is Roj/ 2 in the centered position because (a) the radial velocity distribu- 
tion is assumed to be symmetrical about the midplane, and (b) the circumferential 
velocity is zero at the stator wall and Roo at the rotor wall. Von Pragenau 1 s 
analysis demonstrates that the damper seal yields a lower asymptotic circumferential 
velocity which implies a reduction in the destabilizing cross-coupled stiffness 
coefficient k and a consequential improvement in rotordynamic stability. 

Von Pragenau employs an approximate "short-seal" analysis to develop analytical 
expressions for the rotordynamic coefficients of constant clearance seals. The 
development of these analytical expressions is lengthy and difficult. The combined 
analytical-computational approach used in this development yields an exact 
numerical solution to the governing equations for both constant-clearance and 
convergent-tapered seals with significantly less labor. Following a slight 
modification to Hirs T [7] governing equation to account for different surface- 
roughness conditions on the rotor and stator, the analysis procedure is basically 
that of reference [4] and [5] . 


GOVERNING EQUATIONS 


Figure 2 illustrates a differential element of fluid having dimensions Rd0, 
dz, and H (Z, 0, t) . The upper and lower surfaces of the element correspond to 
the rotor and stator seal elements and have velocities of Ra) and zero, respectively. 
The bulk velocity components of the fluid are Uq and U z ; i.e., these are the averages 
across the fluid film height H of the circumferential and axial fluid velocities. 

The essence of Hirs* formulation is the definition of the wall shear stress x^ as 
the following empirical function of the bulk flow velocity V^ relative to the wall 
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The bulk flow velocities relative to the rotor and stator are, respectively 
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Hence, the shear stress at the rotor and stator are 
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( 4 ) 


Hirs* formulation assumes that the surface roughness is the same on the stator and 
rotor; hence, the same empirical constants mo, no apply to both surfaces* The 
formulation of Eq. (4) accounts for different surface roughnesses in the seal 
elements via the empirical constants (mr, nr), (ms, ns) for the rotor and stator 
surfaces. 


The components of wall shear surface stress in the Z and R0 directions are 
T 
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Summing forces in the Z and R0 directions for the free-body diagram of figure 2 (b) 
yields the following momentum equations*: 
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The bulk-flow continuity equation is 
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*The continuity Eq. (6.c) has been used to simplify these momentum equations. 
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These equations may be nondimens ionalized by introducing the following variables: 


u z = U z /V> u 0 “ V R “’ P = p /pv? 


h = H/C, T — t/T, z - Z/L 


(7) 


T - L/V, b = V/RW 


where C and V are the average clearance and axial velocity, respectively. The 
resultant equations are 

_ ms-KL 


, 3p ns / L 

-7 


R 


ms 


u 


4. BE f L \ mr 

2 u/ • 


1 + 


1 + 


e 


bu 


U 0 -1 


z~\ mr +1 


ibu 


+ h 


3u 3u 

JT + V“ T) TT + 


Z / J 


9u 

U z 


u 


2 u 2 
z 


( 8 ) 


-b 


hf-f (fj K a 


ms 






*• E 


3u fl 3u 0 

3t + u e^ wT ^ air + 


M¥J] 

3u el 

u z 3zJ 


ms +1 

2 


u z u 0 


mr +1 

2-1 2 


u z (u 0- 1) 


8 (bu fl ) 3(hu ) 


PERTURBATION EQUATIONS 


Seal Geometry 

Figure 3 illustrates the geometry for a tapered seal. At the centered position, 
the clearance function is defined by 


H Q (z) = (C + ^) - CtZ = [1 + q (1 - 2 z) ] C = f C 


(9) 


where a is the taper angle, and 
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(10) 


c = (C 0 + C 1 ) / 2 , 



V* 

c 0 + c 


1 

1 


The parameter q is a measure of the degree of taper in a seal and varies from zero, 
for a constant-clearance configuration, to approximately 0.4 for a maximum-stiffness 
seal design [8] . 


Perturbation Analysis 

The governing Eqs. (6) define the bulk-flow velocity components (ug, u z ) and 
the pressure, p, as a function of the spatial variables (R0, z) and time, t. An 
expansion of these equations in the perturbation variables 


u z = u z0 + £ u zl’ h = h 0 + £ h l 

u e = u 9o + e u 0 i> p = p 0 + e Pi 
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where £ = e/C is the eccentricity ratio 


yields the following equations : 


Z eroth-Order Equations: 

(a) Axial-Momentum Equation 
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(b) Circumferential-Momentum Equation 
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(c) Continuity Equation 
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First- Order Equations 


(a) Axial-Momentum Equation 
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(b) Circumferential-Momentum Equation 
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(c) Continuity Equation 
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Most of the parameters of these equations are defined in Appendix A. The quanities 

O , 0 are defined by 
s’ r J 
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(14) 


where the wall friction factors are defined by 
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These expressions correspond to Yamada f s [9] test correlation for flow between 
rotating annulli. 


SOLUTION PROCEDURES 


Zeroth-Order Equations 

The zeroth-order equations define the steady-state leakage and the circumferential 
velocity development Ug Q (z) due to wall shear. The governing equations, Eqs. (12), 
are coupled and nonlinear through the dependency of the coefficients a , a , u 
and V. The equations must be solved iteratively to determine the average lelkage 
velocity V corresponding to a specified pressure drop AP and the circumferential 
velocity distribution u- (z) . The resultant solution defines the leakage coefficient 
C^ of the leakage AP relationship 

AP - C d (16) 

The pressure drop at the entrance is defined by 
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pv 2 (1+g) 
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where \ is an entrance-loss coefficient which is generally on the order of 0.1 to 

0.5. 


320 



First-Order Equations 

The governing first-order equations define Pi( z > T) > u ^ (z, 0, t) , and Uq^ 
(z, 0, t) resulting from the seal clearance functions h^ (0, xf. The clearance 
H is defined in terms of the components of the seal-journal displacement vector 
(X, Y) by 

H = H - X cos0 - Y sin0 (17) 

Hence, by comparison to Eq. (11), 

eh^ = -x cos0 - y sin0 (18) 


where 

x = X/G, y - Y/C 

Note that h^ is not a function of z, and its time dependency arises from the dis- 
placement variables x (t) , y (t) . 

To satisfy circumferential continuity conditions, the following solution format 
is assumed: 
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Substituting from Eq. (19) into Eq. (13) eliminates 0 as an independent variable, 
and yields six real equations. By introducing the complex variables 
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into these equations, the following complex equations are obtained 
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with the A. coefficients defined in Appendix A. The time dependency in these 
equations is eliminated by assuming a harmonic seal motion of the form 

A R 0 jftt jQTx 
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1 C 0 


( 22 ) 


where is a real constant. The associated harmonic solution can then be stated 

A , v — f \ jfiTx 
u zl ^ Z> = u zl ^ e 


u Ql (z, t) = u 01 (z) e^ TT 
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Substitution from Eqs. (22) and (23) into Eq. (21) yields 
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The following three boundary conditions are specified for the solution of 
Eq. (24): 


(a) The exit pressure perturbation is zero; i.e.. 
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P 1 (L) = 0 


(27) 


(b) The entrance circumferential velocity perturbation is zero; i.e., 

u 01 (O) = 0 (28) 

(c) The pressure loss at the seal entrance is defined by 

P s - p(0, 0, T) = | u z 2 (0, 0, T) (1 + 5) 
which yields the following boundary condition: 

P 1 (0) = *(! + £) u zl (0)/(l+q) (29) 

Solution of the differential Eqs. (24) in terms of the boundary conditions is 
relatively straightforward, yielding a solution for the velocity and pressure 
fields of the form 
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Dynamic Coefficient Definitions 


Having obtained the pressure-field solution of Eq. (30), solution for the 
dynamic coefficients will now be undertaken. The reaction-force components acting 
on the rotor due to shaft motion are defined by 
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of Eq. (19) , these integrals further reduce to 
-eRLirpV 2 q / 1 P lc dz; F y (t) = -eRLirp? 2 P lg dz 


(31) 


The motion defined by Eq. (22) is orbital at the precessional frequency Q and 
radius R . This statement may be confirmed by comparing Eq. (18) with Eq. (22) 
to obtain 


X - Cr cosfit, Y = Cr sinQt (32) 

o - o 

Definition of the reaction forces is simplified by performing the integration 
at a time when the rotating displacement vector is pointing along the X axis, i.e., 
when fit = 0. Eq. (23) shows that p^ and p^ coincide for this time and location. 
Hence, Eq. (31) yields the following component force definitions parallel and 
normal to the displacement vector 

F (ST) = -r (TTRLpV 2 ) z 1 f„ (z)dz 
r o o jL 
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F fi (QT) = -r (TTRLpV 2 ) f 1 f, (z)dz 
U o o 3s 

The useful nondimens ional version of these equations is 

F r (OT) 2 L 1 

ttRAPR = * c7~ f 3c ^ z)dz 

o d L 

(33) 

F 0 (ftT) , L 1 

ttRAPR = " cT“ ~jT *o f 3s (z)dz 

o d C 

where R q = Cr is the amplitude of seal motion. The components are expressed as 
function of QT, because, for a given seal geometry (L, R, C) and set of operating 
conditions (AP, 03 ), the excitation frequency £2T is the only independent variable, 
Stated-differently, Eq. (33) provides a frequency-response solution for the reaction 
force components. 

To calculate seal coefficients, a comparable statement of reaction-force 
components is developed from the following nondimens ional statement of Eq, (1) 
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Substitution from Eq. (32) yields 
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Hence, the dynamic seal coefficients (K, k, G, c, M, m) may be obtained by 
comparing the solution obtained by Eq. (33) with Eq. (35). More specifically, 
they are obtained by a least-square curve-fit of the solutions stated on the 
right-hand side of Eq. (35). 


TANGENTIAL VELOCITY DEVELOPMENT 

The frames of figure 4 illustrate the circumferential velocity development 
u 0o (z) which is predicted by Eq. (12b) for the same and different rotor and stator 
surface roughnesses. Roughness is characterized by the empirical coefficients 
(mr, nr), (ms, ns). For the figures illustrated, the roughness of a smooth surface 
is defined by Yamada's [9] coefficients 

mo = -0^25, no = 0.079 

while the parameters 
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mo = -.0024, no = 0.0262, 

corresponding to a theoretical relative roughness £ = e/2C = 0.1 are used for the 
rough surface. Observe that the solution converges towards one half irrespective 
of whether both surfaces are smooth or rough. However, in figure 4 (b) the 
asymptotic solution is less or greater than one half depending on whether one 
uses a smooth-rotor /rough-stator or a rough-rotor/smooth-stator combination. The 
results of figure 4 substantially support von Pragenau's [6] central conclusion 
concerning the desirability of a rough-stator/smooth-rotor combination. 


EXPERIMENTAL RESULTS 


Introduction 

The test results reported here were developed as part of an extended, NASA- 
funded, high-Reynolds-number test program of pump seal configurations in support 
of the SSME (Space Shuttle Main Engine) development program. High-Reynolds numbers, 
which are comparable to those achieved in the cryogenic turbopumps of the SSME, are 
achieved by using CBrF~ as a test fluid. This is a DuPont -manufactured refrigerent 
and fire extinguisher fluid (Halon) which combines a high density and low absolute 
viscosity to achieve very low kinematic viscosity, actually less than liquid 
hydrogen [10]. Details of the flow loop are provided in [11]. 

Figure 5 illustrates the test apparatus. The test fluid enters the center 
and discharges axially across the two test seals. Seal inserts are pressed into 
cylindrical seats in the housing. The rotor segments of the seal are mounted 
eccentrically in the rotor with an eccentricity A. Hence, rotor rotation generates 
a synchronously-precessing pressure field. Axially-spaced, strain-gauge, pressure 
transducers are provided to measure the transient pressure field. Capacitance- 
type proximity probes are provided to simultaneously measure the rotor motion X(t) , 
Y(t) relatives to the housing. The transient pressure measurements are integrated 
to define F^/A and Fq/A, the force coefficients parallel and normal to the seal 
eccentricity vector. In any test, five to ten cycles of data, containing on the 
order of 2,000 data points, are analyzed, yielding a calculated average and standard 
deviation for F /A and F fi /A. The test results reported here were carried out to 
provide answers for the following questions: 

(a) How do predictions from the current theory compare to test results? 

(b) For damper seal configurations, (rough-stator/smooth-rotor) how do 
various roughness treatments compare in terms of leakage, stiffness, 
and damping? 

Stator Configurations 

Tests were carried out on the following stator configurations: 

(a) smooth finish, 

(b) knurled-indentation roughness, 

(c) diamond-grid post pattern. 
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(d) round-hole pattern. 

All seals use the same radial clearance, = .527 ram. Seals b through d are 
illustrated in figures 6. The knurled-indentation roughness pattern is the same 
as that used in current test versions of the HPOTP (High Pressure Oxygen Turbopump) , 
and the seal insert was supplied by Rocketdyne division of Rockwell International, 
the manufacturer of the SSME. The diamond-grid post pattern was manufactured by a 
milling operation which produced grooves which left the square post pattern. The 
round-hole pattern was also produced by a right-hand milling operation. 

Empirical Turbulence Coefficients 

With reference to the adequacy of current analysis, the stator and rotor 
roughness is characterized in terms of empirical coefficients. These coefficients 
must be calculated from the static test data before a theoretical prediction can 
be made for F^/A and Fq/A, and calculation of these coefficients is the subject of 
this subsection. 


In the apparatus of figure 5, a smooth-seal insert is used in the left-hand 
side, while the damper-seal stators were inserted in the right hand side of the 
housing. To the extent possible, the same "very-smooth” finish was provided for 
both the smooth-seal insert and the rotor. Leakage rates and pressure gradients 
were measured for both the smooth and damper seals for all dynamic tests. 


The steady-state axial pressure gradient equation has the form 


3p 




Hence, with a measured pressure gradient and a known density p and axial velocity 
V, the parameter a can be calculated, a is related to the friction-factor co- 
efficient X by 


The smooth-rotor/smooth-stator data were used to calculate O and X values which 
were assumed to apply for both the rotor and smooth stator. r From tfie X versus 
0) and R data, the empirical coefficients mr, nr of the following friction-factor 
formula are calculated 
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(36) 


on a least-square basis, yielding 


nr = 0.0674, 


mr = 


0.217 


For the smooth-rotor/rough-stator combinations, a 
which is related to the corresponding rotor a 


a = 
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a +a 

r s 


combined a is measured, 

Q ' 

and (rough) stator by 


» 2a - 
c 


(37) 


This formula was used to calculate a for the rough stators by using measured values 
for 0 and calculating a value for 0 S from Eq. (36) with the parameters of Eq, (37). 
The empirical coefficients obtained for the stator inserts are provided in table 1, 
together with an estimate of the relative-roughness parameter corresponding to 
pipe-friction data. The results are generally consistent with expectations, 
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except for the positive value for ms obtained for the hole-pattern stator; however, 
over the Reynolds-number range tested, the combined friction-factor actually 
increased (slightly) with increasing R^for this stator insert. 


Dynamic Test Data 

For a given seal configuration, a test matrix is obtained by varying the axial 
Reynolds number and running speed. The R range varies between the maximum flow 
capacity of the supply pump and the minimum AP sufficient to generate reasonable 
transient pressure signal amplitudes. For a given R & value, the running speed is 
varied sequentially over the running-speed capacity o? the drive motors. Figures 
7 through 10 illustrate theoretical and experimental results for the four stators 
tested. An inspection of these results demonstrates "reasonable" agreement between 
theory and experiment for Fq/A but much larger F r /A magnitudes at lower speeds than 
predicted. Further, the magnitude of F r '/A decreases more rapidly with increasing 
running speed than predicted. 


DISCUSSION OF EXPERIMENTAL RESULTS 


Comparison to Theory 

If a circular orbit of the form 

X = A cosajt Y = A sinoot 

is assumed, Eq. (1) yields the following definition of force coefficients 

F /A = -K - COJ + Mu) 2 
r 

Fg/A - k - Co) 

where the cross-coupled mass coefficient m has been dropped as being negligible in 
comparison to the influence of k and C. At first glance, these equations suggest 
that sufficient independent equations could be obtained, in the present apparatus, 
to independently calculate all the rotordynamic coefficients by holding the flowrate 
constant and running at three different speeds. However, the. fact that the co- 
efficients depend on a) precludes this approach. While K, C, and M are weak functions 
of a) through their dependence on a, the "cross-coupled" coefficients k and c are 
linear functions of a). In fact, if the fluid is prerotated prior to entering the 
seal such that the inlet tangential velocity is Uq q (o) = Ro)/2, then theory predicts 
that k = Co)/2, c * Ma), and 

F /A = -K, F q /A = -Co)/2 
r y 

The present test apparatus provides no intentional prerotation, and the expected 
result is of the form 

k * b^Cto/2, b^<l 

c = b 2 Ma), b 2 <l 

Fq/A— -C ef fcd = -C(l-b 1 /2)0) 
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F r /A 3£ -K ef + H ef C0 2 - -K + M(l-b 2 )0) Z 

The term denotes the "net damping coefficient" resulting from the drag force 
Co zk and the forward whirl excitation force kA. A direct comparison between theory 
and experiment is obtained by curvefitting the theoretical and experimental results 
for the F^/A and Fg/A to obtain predictions for K^, C and * Note that the 
procedure of curvefitting the data with respect to U) e?Iminates e the running-speed 
dependency. Further, K ^ is the zero-running speed intercept of the F r /A versus oo 
curve, and is the slope of the Fq/A versus to curve. r 

A comparison of measured and experimentally-derived values for K^, C^, and 
are given in table 2 for the stators tested, and support the following general 
conclusions: 

(a) Direct stiffness values are substantially underpredicted by theory. This 
result is consistent with earlier water test results [5, 12]. Improved 
correlation generally results at larger C/R ratios. 

■(b) Net damping coefficients are overestimated by theory, but the agreement 

is reasonable and generally improves with increasing R . . 

ao 

(c) The added-mass coefficient is substantially underpredicted by theory. 
However, this result is at variance with earlier water-test results 
[12] which show an overestimation of the added-mass coefficient. 

Relative Performance of Stators 


Dynamic Coefficients 

Figures 11 and 12 illustrate K f and C ^ for the stators tested versus AP, and 
can be used for direct comparison of the stiffness and effective damping of the 
roughness designs. The results support the following conclusions. 

(a) The knurled-indentation and the diamond-grid stators are, respectively, 
the most and least stiff. The hole-pattern and smooth stators have 
comparable stiffness. 

(b) The hole-pattern and diamond-grid stators provide, respectively, the most 
and least net damping. The smooth and knurled-indentation stators have 
comparable net damping coefficients. 


The disappointing performance of the diamond-grid stator is related to its larger 
average clearance. The relieving operation which yields the posts yields an average 
clearance of 0.889 mm as compared to the 0.527 mm minimum clearances of the remaining 
configurations . 


Leakage Performance 


To evaluate leakage performance, 
coefficient definition 


AP = C d 


pF 

2 


C T is defined using the conventional discharge 
Li 


which yields 
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Hence , 



The coefficient is a nondimens ional relative measure of the leakage to be 
expected through seals having the same radius. Figure 13 illustrates C versus 
AP for the seal stators and demonstrates that the round-hole pattern ana smooth 
stators have, respectively, the best and worst performance. The knurled- indentation 
pattern has a slightly better leakage performance than the diamond-grid pattern. 


CLOSURE 

A theory is presented, based on a simple modification of Hirs 1 turbulent 
lubrication equations, to account for different but direct ional ly-homo geneous 
surface roughness treatments for the rotor and stator of annular seals. The 
theoretical results agree with von PragenauVs predictions that a "damper seal" 
which uses a smooth rotor and a rough stator yields more net damping than a 
conventional seal which has the same roughness for both the rotor and stator. 

Experimental results for four stators confirm that properly-designed 
roughened stators can yield higher net damping values and substantially less 
leakage than seals with smooth surfaces. The best seal from both damping and 
leakage viewpoints uses a round-hole-pattern stator. Initial results for this 
stator suggest that, within limits, seals can be designed to yield specified 
ratios of stiffness to damping. Additional testing for this type of seal is 
scheduled for 1984-1985 to examine the influence of hole depth, hole shape, and 
the ratio of hole-relieved area to total surface area. 


APPENDIX A: PERTURBATION COEFFICIENTS 
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TABLE I. - EMPIRICAL TURBULENCE COEFFICIENTS MS, NS, 
AND ESTIMATES FOR RELATIVE ROUGHNESS 



ms 

. . 

ns 

e = e/2C 

Smooth 

-.240 

.0989 

.00069 

Knurled-Indentation 

-.136 

.0697 

.022 

Diamond-Grid 

-.0350 

.118 

.460 

Hole Pattern 

.0190 

.0150 

.058 
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TABLE II. - A COMPARISON OF THEORETICAL AND MEASURED VALUES FOR 
EFFECTIVE STIFFNESS, DAMPING AND ADDED-MASS COEFFICIENTS 


SEAL 

R ao 

■H 

MM 

MEF 

ex 

MEF 

th 


120,100 

2.182 

0.6636 

1 . 38 

Smooth 

170,600 

2.216 

0.6900 

1 .38 


270,500 

1 .297 

0.7567 

4.45 


385,200 

1.142 

0.9350 

6.86 


495 ,700 

1.261 

0.9964 

9.35 


500,900 

1 .478 

0.9401 

0.986 



1 .630 

0.5613 

mm 

Knurled- 

160,200 

1 .598 

0.8166 

Efffl 

Indentation 

335,600 

1.678 

0.9836 

18,6 


Esifiiisfl 

1.664 

1 .027 

1 1 .2 


368,200 

1 .706 

1.041 

4.23 



1.745 

0.4265 

BP 

Diamond Grid 

159,800 

1.310 

0.4795 

2.61 


iBH 

1 .065 


7.58 


349,900 

0.8026 

0.9583 

1 1 .9 



1.128 

1.012 

4.83 


89 , ! 41 0 

1 .573 

0.4506 

1 .94 

Hole Pattern 

130,400 

1 .768 

0.6997 

2 .25 


159,700 

1 .876 

0.7867 

2.58 



2.175 

1 .205 

■ns 
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Figure 1. Neck-ring and interstage 
seals. 


Figure 3. Tapered seal 
geometry. 




SEJ& 

39 2 


)HdZ 



Figure 2, Seal differential element. 



P ROUGH STRTOR/ROUGH ROTOR 
O SMOOTH STRTOR/SMOOTH ROTOR 
Ra- 2SI4R.5 Rc-207254. 



fix l is I Location 

Figure 4(a). Predicted circumferential velocity 

development for the same rotor and stator roughnesses. 


□ SMOOTH STRTOR/ROUGH ROTOR 
O ROUGH STRTOR/SMOOTH ROTOR 
Ra- 2S14R.S RC-2072S4, 



fix t a I Location 

Figure 4(b), Predicted circumferential velocity 

development for different rotor and stator roughnesses. 
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FCTANGENTIA1DVA (HN/'M) F ( RAO I AL ) /A C MN/M ) 



Figure 6(b). Diamond grid, 
post-roughness stator 
insert . 



Figure 6(c). Round-hole 
pattern stator insert. 


FINITE-LENGTH THEORY 
SMOOTH STATOR / SMOOTH ROTOR 
C( RADIAL) - .5271 mm 

Mapa-jrael » ■ Thacry 



a Ro - 8A860.2 
O Ra - 120085. 
A Ro - 170565 . 
+ Ro • 270556* 
O Ro - 366172 » 
H Ra - 4S5726 , 
n Ra • 500666. 


finite-length theory 

SMOOTH STATOR ✓ SMOOTH ROTOR 
Ct RAO I A l) • .6271 mm 



□ Ra - 6A060.2 
O Ra - 120265. 
A Ra - 170&8S. 
+ Ra - 270556. 
O Ra - 385172, 
H Ro - 495725, 
U Ra - 600966. 


Figure 7. 

Measured and theoreti 
results for F /A and 
smooth stator* 
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FI N I TE- LENGTH THEORY 
ROCKETOYNE STATOR / SMOOTH ROTOR 
C ( RADIAL ) - .6271 mra 


\ 

Z — — — Measured " 1 Theory 



ROTOR SPEEO (RAO/SEC) 

S Cl RADIAL.) - .5271 mm 

z 

w — — ■— h©o*ured — — Theory 



D Ra - 6S860 , 2 
O Ra - I 1 5070 . 
A Ra - 160198. 
+ Ra “ 335625. 
O Ra - 3S04A7. 
K Ra - 388227. 


Q Ra - 6A860.2 
0 Ra • 116670. 
A Ra - 160188, 
+ Ra - 335625. 
O Ra - 3S04A7 , 
H Ra - 369227. 


Figure 8. 

Measured and theoretical results 
for F /A and Fq/A; knurled- 
indenlation stator. 
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FINITE- LENGTH THEORY 
DIAMOND GRID STATOR / SMOOTH ROTOR 
C< RADIAL -AVERAGE) - .88 A mm 


M««*ur«d Thaory 



ROTOR SPEED ( RAD/SEC) 

C( RADIAL- AVERAGE \ - ,88A mm 



ROTOR SPEED (RAO/SEC) 


□ Ra • 7052 A. 
O Ra - US087, 
A Ra - 15A848. 
+ Ra - 355625 . 

0 Ra - 34AA00. 

1 Ra • 383AS0 . 


D Ra - 7852A. 
O Ra - l 1 5087 . 
A Ra - 1SA848. 
+ Ra - 355625, 
O Ra - 34AA00. 
B Ra - 385A50 . 


Figure 9. 

Measured and theoretical 
results for F^/A and F^/A; 
diamond-grid stator. 
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F t TON GENT 1 R*L ) / f\ (MN/M) 


FINITE-LENGTH THEORY 
ROUGH STATOR i HOLE PATTERN 1 


M«e*ur«d Theory 



ROTOR SPEED ( RAO/SEC ) 
ROUGH STATOR i HOLE PATTERN 1 


0 Ro • 0A4 14*8 
O Ro • 130385. 
A Ro - 1SA727. 
+ Ro - 320438. 


Figure 10, 

Measured and theoretical 
results for F /A and Fg/A 
round-hole pattern. 


Measured 


Theory 



ROTOR SPEED ( RAD/SEC ) 


p 

R« - 

8R41 4 , 6 

o 

Ro - 

1 30382 , 

A 

Ro • 

15A727. 

+ 

Ro - 

328438, 


P KNURLED 
O SMOOTH 
A DIAMOND 
+ ROUND - 


INDENTATION (C- .5271 mm) 

(C- .5271 mm) 
GRID (Cm- . 88A mm) 

HOLE 



Figure 11. K gf versus AP for stator inserts. 
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